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Abstract 



£T) ■ We investigate a scenario where conditions for Mott dissociation of 

CN , charmonium are fulfilled locally in ultrarelativistic heavy ion collisions 

at the short time scales before hadronization and study possible conse- 
quences to be drawn from observation of "anomalous" thresholds in the 
J/ip suppression pattern. We suggest that for a compilation of exper- 
iments with different kinematics, variables presently in use to display 

C\| ' the observed suppression pattern should be rescaled by the 7 factor of 

the cc state relative to the participant center of mass. This systematics 

qv ', can be tested by inverse kinematics experiments as, e.g., Pb-S collisions 

and by sampling data sets in px and y of the dimuon pairs. 

"o ! 1 Introduction 

Recently, an "anomalous" threshold behaviour in the centrality dependence 
of J /if) suppression has been observed by the NA50 experiment [U in Pb- 
^ Pb collisions at 158 A GeV/c. It has been speculated that this indication 

of a critical behaviour might signal quark-gluon plasma formation 0, || as 
suggested earlier by Matsui and Satz |] to be a consequence of the Mott 
mechanism for charmonium bound states in a plasma. 

Properties of bound states are modified in dense matter. Above critical 
densities timott (temperatures Tmott) they undergo a Mott transition to un- 
bound (but correlated) states in the continuum. For ground state hadrons 
as light quark bound states this transition corresponds to the phase transition 
from hadronic to quark matter || . For the deeply bound J /if) at this transition 
there is still an energetic threshold to be overcome either by kinetic hadron or 
parton impact or by further compressing and heating the plasma [|, |7|, ||, ||. 
In this contribution, we want to investigate a scenario where conditions for 
Mott dissociation of charmonium are fulfilled locally at the short time scales 
before hadronization sets in and study possible consequences for the NA50 
experiment. 



2 Mott dissociation of heavy quarkonia [] 

In order to determine the critical density riMoTi{QQ*) at which the QQ* pair 
will not evolve into quarkonium but rather dissociate into a pair of heavy 
mesons, we employ here the heuristic criterion that at this density the nearest 
neighbor for a given heavy quark Q is not a heavy antiquark Q but a light 
antiquark q of the medium H . According to elaborated quantum mechanical 



approaches [11, [L2|, at short time scales r < n after the creation and color 



neutralisation, the heavy quark pair evolves like a diverging wave packet 

*QQ 



< r 2 > QQ , (r) =< r 2 > qQ (r/r^, (1) 



until it reaches the asymptotic value < r 2 >qq of a quarkonium state in free 
space, we take f3 — 2 |9j, fL3[l . The mean squared distance of a heavy-light quark 



pair is obtained from the nearest neighbor distribution function |7| for a plasma 
density which is proportional to the participant density n part (t) measured in 
the c.m.s. proper time 

< r 2 > Q g* (r) = const /n paTt (rj) , (2) 

where r is the time in the rest system of the quarkonium* and the 7 factor 



j =^cosh 2 (Ay)+ p 2 T /M 2 + ^ (3) 

takes into account that the QQ* state is Lorentz boosted relative to the c.m.s. 
of the participants. The difference in the 7- factors of S-U and Pb-Pb collisions 
to be discussed in Sect. 3 is mainly due to different central rapidities of about 
y c = 2.4 and y c = 3.0, respectively, in the rapidity shift Ay = y c — y fJ ,+ IJl -. 

During the phase when n part (t) rises to its maximum at t = t q j ~ 1 
fm/c, it eventually exceeds the critical Mott density umott^QQ*) for which 
< r 2 >Qq*=< r 2 >qq* and instead of a quarkonium state a pair of heavy 
mesons will emerge from that region in the reaction plane for which 

n part (b,s) > nuoTr(QQ*) = 7 /3 ™mott(QQ) , (4) 

where n MoTT (QQ) = const (rf/to) 13 / < r 2 >qq is the critical Mott density for 
a QQ* pair created at rest in the medium (7 = 1). In Fig. 1, we demonstrate 
for /3 = 2 that the occurence of "anomalous" behaviour (Mott effect) in a 
given nucleus-nucleus collision depends on the kinematical conditions. The 
dependence on the scaled variable n paiTt /'y 13 , however, is expected to uncover 
the threshold behaviour due to the Mott dissociation. 



*The asterix indicates that the QQ state is a color singlet but not yet a state of the 
charmonium spectrum and will be denoted as in statu nascendi, see jlOf . 
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Fig. 1. Criterion for Mott-dissociation of charmonium states. Central par- 
ticipant densities (long dashed lines) as a function of the impact parameter b 
(bottom scale) or equivalently the transverse energy Et (top scale) have to be 
scaled by the 7- factor of the charmonium Lorentz-boost relative to the c.m.s. 
of the collision, since the sizes of the charmonia in statu nascendi at tq depend 
on 7. Solid lines correspond to the present data of the S-U (left panel) and 
the Pb-Pb (right panel) experiments, resp. The horizontal dashed lines show 
the critical Mott densities which mark the onset of "anomalous" suppression 
as observed in the Pb-Pb experiment. 



3 Glauber model for anomalous J/ifj suppression 

Previous analyses of J/ip production in pA as well as AB collisions within the 
Glauber model have shown [||, [14]] that nuclear absorption gives a satisfactory 
description of J/-0 suppression data except for the recent Pb-Pb data sample 
[|I|] which therefore has been termed "anomalous". The relevant quantity in 
this description is the survival probability S'nuciW = S d 2 s S nuc \{h, s) for a J/ip 
formed in an AB collision at impact parameter b — |b| which is an integral 



over the survival probability density Wt 

q (u ,_ (1 ~ (1 ~ a ahs T A (s)) A )(l - (1 - a ahs T B (s - h)) B ) 
nUCl( ' j " <*, AB f cPs T A (s)T B (s - b) {b) 

in the transverse plane. T A (T B ) is the nuclear thickness function for a Woods- 
Saxon density distribution of the target (projectile) nucleus normalized to 
unity. A cross section cr abs = 7.3 mb which is compatible with pA data gives 
an excellent description of J/ip production in AB collisions as performed by 
the NA38 experiment. For Pb-Pb collisions the nuclear absorption (^) alone 
fails to describe the data, see Fig. 2. Additional comover absorption |L5| does 



not give an overall satisfactory description of the whole data sample even if a 
continuous increase of the absorption cross section with n paTt is assumed, see 
also ||. 

The rather pronounced threshold observed in the Pb-Pb data hints at a 
critical behaviour as, e.g., predicted by the Mott dissociation effect con- 
sidered in the previous Section. We want to implement this effect into the 
Glauber model approach by cutting those regions off the reaction plane where 
the Mott condition (01) is fulfilled. Considering a 30% feeding of the observed 
J /ip from decays of Xc and neglecting the small correction for ip' we obtain the 
total J/ijj survival probability 

S(b) = 0.7 S 4 ,{b) + 0.3 S x (b) 

S t (b) = /d 2 s5 nucl (b,s)eK art (b,s)- 7 /3 n M0TT (z)]. (6) 

According to the Glauber model the density of participants in the trans- 
verse plane for a given impact parameter is n part (b,s) = AT A (s)(l — (1 — 
<T in eiT B (s - h)) B ) + BT B (s - b)(l - (1 - a- mcl T A (s)) A ) . The total survival 
probability @ is a generalization of equations employed in Refs. 0, |3|] . The 
main difference is that in Eq. @ we account for the fact that suppression of 
quarkonia in statu nascendi depends on the time scales involved. Consequently, 
the critical Mott-densities (|4|) have to be rescaled with the 7-factor (|3|) which 
implies a dependence on the rapidity shift Ay and the transverse momentum 
Pt at which the dimuon pair is detected. In order to use different experimental 
situations for a more precise determination and identification of the thresholds, 
it is necessary to define a variable respective to which the threshold behaviour 
appears universal. After inspection of Eq. (P) we propose here the scaled vari- 
able < n part (b) > /7 /3 , where the mean number of participants is defined [[| by 
< ?Vrt(fr) >= / d2s n hard(b, s) n part (b, s)/ / cPs n hard (b, s) , and the density of 
hard collisions is n haiTd (b, s) = A B T A {s)T B {s — b)a^ N . The numerical results 
are shown in Fig. 2. 
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Fig. 2. J/ip suppression in nucleus-nucleus collisions normalized to pure nucle- 
onic absorption. Whithin our approach, all data should approximately follow 
a universal curve as function of the averaged participant density divided by 
7^, where the two thresholds correspond to the onset of Mott dissociation 
for the Xc and the tp state, respectively. The open diamonds show an estimate 
for the result of an inverse kinematics experiment Pb-S which fills the gap 
between present S-U (filled squares) and Pb-Pb (filled circles) data sets. 



4 Conclusions 

We have reanalyzed recent data on J/ip production in S-U and Pb-Pb colli- 
sions at CERN-SpS and found that besides conventional sources an additional 
absorption mechanism with critical behaviour is necessary. In the present work 
we have generalized the Glauber model approach to J/ip suppression by imple- 
menting the critical Mott dissociation effect for charmonia in statu nascendi. 
There are new aspects in this analysis: (1) there are two thresholds ex- 
pected to occur in the "anomalous" J/ip suppression pattern due to the criti- 
cal Mott effect for the lp (% c ) and the Is (J/VO state; (2) the corresponding 
threshold values for impact parameter or transverse energy depend on the kine- 
matics of the cc pair relative to the hot and dense medium which it experiences 
in statu nascendi; (3) conclusions from J/ip suppression for the diagnostics of 
the state of matter (QGP formation) depend on the kinematics. As a conse- 
quence, QGP might be present in S-U collisions although J/ip production as 



measured by NA38 at < 7 2 >= 2.6 is not "anomalous". 

In order to investigate the threshold behaviour of J/ip suppression we have 
proposed the variable < n paTt > /j", where (3 = 2 is supported by present 
data. We find it most useful to perform systematic experimental studies by 
(1) using inverse kinematics (for lead beam on sulphur, see estimate in Fig. 2) 
and (2) sampling low vs. high p? and/or low vs. high rapidity bins. These 
analyses can be performed using data as provided by the NA50 experiment. 

The theoretical approach to the kinetics of the hadronization process for a 
heavy quark pair in a dense plasma has been simplified here by the introduction 



of an instantaneous breakup. A more detailed investigation is in progress ||T7 |. 
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